A measurement by BOOMERANG of multiple peaks in the 
angular power spectrum of the cosmic microwave background 
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ABSTRACT 

This paper presents a measurement of the angular power spectrum of the Cosmic Microwave 
Background from I = 75 to I = 1025 (~ 10' to 2.4°) from a combined analysis of four 150 GHz 
channels in the BOOMERANG experiment. The spectrum contains multiple peaks and minima, 
as predicted by standard adiabatic-inflationary models in which the primordial plasma undergoes 
acoustic oscillations. These results, in concert with other types of cosmological measurements 
and theoretical models, significantly constrain the values of f^tot, Qbh 2 , Q c h 2 and n s . 



Subject headings: Cosmic Microwave Background Anisotropy, Cosmology 



1. Introduction 

The presence of a harmonic series of "acoustic" 
peaks in the angular power spectrum of the cos- 
mic microwave background (CMB) was predicted 
as early as 1970 (Sunyaev & Zeldovich 1970; Pee- 
bles & Yu 1970). These peaks arise due to the evo- 
lution of pressure waves in the pre-recombination 
universe and are a generic feature (Bond and Ef- 
stathiou 1987) of most, but not all, cosmological 



models (e.g., Hu et al. (1997)). Specifically, a well- 
defined set of peaks is predicted for both adiabatic 
and some classes of isocurvature models of struc- 
ture formation, but not for models that rely on 
topological defects. 

Since the COBE measurement of the ampli- 
tude of fluctuations in the cosmic microwave back- 
ground at the largest scales (Bennet et al. 1996), 
a large literature has developed which predicts, in 
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the context of adiabatic cold dark matter (CDM) 
models, the relative position and amplitude of 
these peaks for different values of the fundamental 
cosmological parameters. A general prediction of 
these models is the presence of a dominant funda- 
mental peak at an angular scale w 1° (£ w 200), 
decreasing in angular scale when f2 decreases. 
Data from a variety of experiments (Miller et al. 
1999; Mauskopf et al. 2000; Hanany et al. 2000) 
including a small fraction of the data from the 
BOOMERANG 1998/1999 Long Duration Bal- 
looning (BOOM/LDB) campaign (de Bernardis et 
al. (2000) ; BOO hereafter) clearly show this feature 
and provide strong evidence for a low curvature 
universe, a generic prediction of many inflation 
models. 

There is also convincing evidence that the 
broad-band average of power at smaller angular 
scales gradually declines in a manner consistent 
with adiabatic CDM models (Padin et al. 2001; 
Church et al. 1997). However, these experiments 
do not have the necessary combination of sensi- 
tivity and sky coverage to convincingly detect or 
reject the presence of harmonics of the fundamen- 
tal peak in the power spectrum. The detection 
of such harmonic peaks would provide strong ev- 
idence that the scenario of structure formation 
from acoustic oscillations in the primordial plasma 
is accurate. 

Present here is an analysis of a larger set of data 
than previously released from the BOOM/LDB 
(de Bernardis et al. 2000) experiment which shows 
clear evidence of multiple peaks in the angular 
power spectrum of the CMB. Data from four sep- 
arate detectors that each observe 1.8% of the sky 
are combined. A new data analysis algorithm 
is used and refined estimates of the beam shape 
and overall experimental calibration are presented. 
The spectrum is consistent with low spatial cur- 
vature, A-dominated adiabatic CDM models. 

2. Instrument 

BOOMERANG is a Long Duration Balloon 
(LDB) experiment designed to measure the an- 
gular power spectrum of the CMB at degree and 
sub-degree scales. For a complete description of 
the instrument see Crill et al. (2001) and Piacen- 
tini et al. (2002). Instrument characteristics are 
summarized in Table 1. 



BOOMERANG was launched for its first LDB 
flight on December 29, 1998 from McMurdo sta- 
tion, Antarctica, and acquired 257 hours of data 
from an altitude of w 39km. BOOMERANG 
is comprised of a 1.2m off axis parabolic mirror 
which feeds a cryogenic mm-wave bolometric re- 
ceiver. Observations are made simultaneously in 
four unpolarized bands centered at 90 GHz, 150 
GHz, 240 GHz and 410GHz. 

The telescope is steerable in azimuth by mov- 
ing the entire gondola, and in elevation by mov- 
ing an inner frame containing both the receiver 
and the optics. The illumination of the optics is 
not modulated by the scan, which minimizes scan 
synchronous optically generated offsets. Exten- 
sive shielding permits observations in the azimuth 
range ±60° from the anti-sun direction, for all sun 
elevations experienced in the antarctic ballooning 
environment. 

3. Observations and Calibration 

Observations arc made by scanning the tele- 
scope in azimuth by 60° peak-to-peak at an an- 
gular velocity of 2°/s (hereafter 2dps) (for the 
first half of the flight) or l°/s (hereafter ldps) 
(for the second half of the flight) at fixed eleva- 
tion. Each day, the elevation is shifted. Observa- 
tions are made at elevations of 40°, 45°, and 50°. 
The scans are centered well away from the Galac- 
tic plane. Interspersed in the CMB observations 
are observations of selected point sources near the 
Galactic plane. The CMB sky coverage is shown 
in Figure 1. 

The results in this paper contain data from the 
entire BOOM/LDB flight (as compared to the sec- 
ond half only in BOO) and incorporate data from 
4 150 GHz channels (as opposed to 1 for BOO). 
This was enabled by using the new data analysis 
techniques described in Section 4 and in Hivon et 
al. (2001). In addition, the most significant source 
of systematic uncertainty at multipoles, I > 400 
is the effective beam size (dominated by pointing 
uncertainties). Since the release of BOO, the under- 
standing of the beam and pointing has been sig- 
nificantly improved, allowing computation of the 
power spectrum out to higher multipoles. 

Observations of extra Galactic sources in the 
main map and scanned observations of bright HII 
regions near the Galactic plane are used to esti- 
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Table 1 
Instrument Characteristics. 



Channel 


Band (GHz) 


NETcmb QiKy/a) 


Beam FWHM (') 


B150A 


148.0 - 171.4 


130 


9.2 ± 0.5 


B150B 


145.8 - 168.6 


Variable 


9.2 ± 0.5 


B150A1 


145.5 - 167.3 


231 


9.7 ± 0.5 


B150A2 


144.0 - 167.2 


158 


9.4 ± 0.5 


B150B1 


144.2 - 165.9 


196 


9.9 ± 0.5 


B150B2 


143.7 - 164.3 


184 


9.5 ± 0.5 


90 (2 Chs) 


79 - 95 


140 


18 ± 1 


240 (3 chs) 


228 - 266 


200 


14.1 ± 1 


410 (4 chs) 


400 - 419 


~ 2700 


12.1 ± 1 



Note. — Summary of relevant instrument characteristics. Only results from 
the 150GHz channels arc presented in this paper. B150B is not used due to 
non-stationary detector noise. The bandwidth limits arc computed to include 
68% of the total detected power for a flat spectrum source. The NET is 
computed at 1Hz. 
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Fig. 1. — Sky coverage. The upper panel shows 
the BOOMERANG 150GHz map. The locations 
of the three bright quasars are circled. The sky 
subset used in BOO (rectangle) and this paper (el- 
lipse) are shown. The bottom panel shows the 
integration time/pixel. The ellipse for this anal- 
ysis was chosen to include the well sampled sky, 
and to avoid the poorly sampled sky 



mate the beam full width half maxima (FWHM) 
for each channel. These values are presented in 
Table 1. While these observations provide good 
statistics on the width of the main lobe, the sig- 
nal to noise of these data are not sufficient to fully 
characterize the near sidelobe response and devi- 
ation from gaussianity of the beam. To obtain a 
detailed model of the near sidelobe response of the 
telescope, a physical optics calculation of the beam 
shape is performed for each channel using the mea- 
sured position of the horns in the focal plane. To 
check the precision of the model calculations, a 
comparison with near field beam maps is made. 
Azimuthal cuts through the telescope beams are 
measured to the ~ —20 dB level. While the 
channels utilizing conical feed horns compare well 
with the beam map data, a ~ 10% discrepancy 
is observed between the calculated and measured 
FWHM of the photometer beams. In all cases, 
the measured FWHM are larger than the calcu- 
lated beam size. This discrepancy is attributed to 
the multi-moded nature of the Winston concen- 
trators utilized by these channels. For these chan- 
nels the FWHM are scaled to fit the measurements 
of RCW38. To correct for the extended nature of 
the source, the angular extent of RCW38 was mea- 
sured independently by the ACBAR instrument 
to be 2.5' (Acbar 2001). The two dimensional far- 
field radiation patterns from the physical optics 
calculation are then used to generate the window 
functions for each channel. 
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The telescope pointing is determined using a 
combination of rate-gyros and an azimuth sun sen- 
sor. In order to recover the long time scale pendu- 
lations, the gyroscopes are integrated with a 400s 
time constant. Given the quoted noise in the gy- 
roscopes of less than 8' /V hour, a pointing uncer- 
tainty of less than 2.7' (la) is predicted. The 
Galactic plane point source observations give a 
pointing uncertainty of 2.5' (la). In the analysis, 
the calculated beam is convolved with the Gaus- 
sian approximation of the pointing uncertainty. 

The pointing solution has been improved in this 
analysis (compared with BOO) by better use of the 
pitch and roll rate gyros, and the regression of a 
thermally dependent offset in the elevation. A re- 
analysis of the pointing jitter in the BOO pointing 
solution, utilizing the apparent centroids of point 
sources along the Galactic plane yields an effective 
beam size of 12.9 ± 1.4', for the original pointing 
solution (as compared to the quoted 10' ± 1' used 
in BOO). When the calibration and beam uncer- 
tainties are taken into account, the new results lie, 
bin by bin, within the overall uncertainty of the 
BOO spectrum, which was restricted to i \/= 600. 
The BOO results are systematically lower at high 
I than those presented here, due to the smaller 
effective beam that was assumed for the BOO anal- 
ysis. However, it is reassuring to note that cor- 
recting the BOO spectrum with the new estimate 
of the BOO pointing jitter eliminates any residual 
discrepancy between the new results and the BOO 
spectrum (de Bernardis ct al. 2002). 

The gain calibrations of the 150 GHz channels 
are obtained from observations of the CMB dipole. 
The data are high-pass filtered with a filter de- 
scribed by F = 0.5 (l - cos for < v < f 
and 1 elsewhere. In order to retain more large- 
scale information than is needed in the anisotropy 
analysis, / is set to the relatively low value of 
0.01 Hz. To compare with the data, we artifi- 
cially sample the CMB dipole signal (Lineweaver 
ct al. 1996), corrected for the Earth's velocity 
around the sun (Stumpff 1980), according to the 
BOOMERANG scanning, and filter this fake time 
stream in the same way as the data. The ldps data 
is then fit simultaneously to this filtered dipole, a 
similarly filtered dust emission model (Finkbcincr 
ct al. 1999), an offset and the BOOMERANG 410 
GHz data for all data more than 20° below the 



Galactic plane. The dipole calibration numbers 
obtained with this fit are robust to changes in 
Galactic cut, and to whether or not a dust model 
is included in the fit; this indicates that dust is not 
a serious problem for the contamination. They are 
insensitive to the inclusion of a 410 GHz channel in 
the fit, which is a general indication that there is 
no problem with a wide range of systematics such 
as atmospheric contamination, as these would be 
traced by the 410 GHz data. 

Overall, the calibration of the spectrum has 
risen by 10% in Q (5% in AT) compared to BOO 
due to a refinement of the dipole calibration (im- 
proved treatment of the time stream filters) and is 
further raised by a better calculation of the beam 
sidclobes. 

4. Data Analysis 

The data are analyzed in four basic steps: i) 
the reduction of the raw data into pointed and 
flagged time streams, ii) the estimation of the 
noise via an iterative map-making algorithm, iii) 
the estimation of the angular power spectrum via 
Monte-Carlo calibrated spherical harmonic trans- 
forms, and iv) estimation of parameters by like- 
lihood methods. The size and nature of the 
BOOMERANG data have required the develop- 
ment of new techniques. 

In the reduction of the raw bolometer data, the 
filter response of the detector and associated elec- 
tronics is deconvolved from the time stream, and 
transient phenomena (predominantly cosmic rays) 
are flagged and replaced in the time stream with 
a constrained realization of the noise. Similarly, 
the RA/Dec pointing for each channel is recon- 
structed from the rate gyros, azimuth sun sensor 
angle, GPS co-ordinates of the package, and the 
focal plane geometry. Details can be found in Crill 
ct al. (2001). 

The receiver noise for each channel is estimated 
from the raw time stream by iteratively solving 
simultaneously for the detector noise spectrum 
n(f) and the maximum likelihood CMB map, 
A = (P t N" 1 P)- 1 P t N" 1 d. The algorithm used 
is an approximate Jacobi method: 
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loop on j 

• nO) = d PA (j) => NO)- 1 = (nnt)- 1 

• A(j +1 ) - £0) = (ptp) _1 ptN(j)- 1 n(j) 
end loop 

where is N is the time-time correlation matrix, 
P is the pointing matrix, d is the time-stream 
data, and n is the noise time stream. If the noise is 
stationary, then N _1 is diagonal in Fourier space, 
and multiplication by N -1 is just a convolution. 
And multiplication by (P^P) P^ represents bin- 
ning into pixels and dividing by hits per pixel. For 
details see Prunet et al. (2000). A complete map 
and noise spectra takes about 15 minutes on an 
alpha-cv67 computer. 

Two closely related estimators were used to 
recover the underlying CMB power spectrum 
Ct from the data. Both methods are based on 
the Monte Carlo Spherical Harmonic Transform 
(MASTER) technique described in Hivon ct al. 
(2001). MASTER allows fast and accurate deter- 
mination of Ct without performing the extremely 
time consuming matrix-matrix manipulations that 
characterize exact methods and limit their appli- 
cability (Borrill (1999)). It can be summarized as 
follows. The spherical harmonic transform of a 
naively binned map of the sky is calculated using 

1/2 

a fast 0(N p - x I) method based on the Healpix tes- 
sellation of the sphere (Gorski et al. 1998). The 
angular power in the noisy maps, Ct, can be re- 
lated to the true angular power spectrum on the 
full sky, Ct, by the effect of finite sky coverage 
(Mtt>), time and spatial filtering of the maps (Ft), 
the finite beam size of the instrument (Bt), and 
instrument noise (Nt) as 

(Ct) = M u>Ft>Bl (Ct') + M "' F ^ ( N *) • 
t' v 

(1) 

The coupling matrix M U ' is computed analyt- 
ically, Bt is determined by the measured beam, 
Ft is determined from Monte-Carlo simulations of 
signal-only time streams, and Nt from noise-only 
simulations of the time streams. 

The simulated time streams are created us- 
ing the actual flight pointing and transient flag- 
ging. The signal component of the simulated time 



streams is generated from simulated CMB maps 
and the noise component from realizations of the 
measured detector noise n(f). Ft and Ng are de- 
termined by averaging over 150 and 200 realiza- 
tions respectively. Once all of the components are 
known the estimation is carried out in two ways. 

In the first case the power is determined by solv- 
ing directly for the unbiased estimator Ct of eqn. 1. 
The uncertainties in the estimated top-hat binned 
Ct spectrum are measured by averaging over en- 
sembles (typically 400 realizations) of signal+ noise 
simulations created using a best fit model power 
spectrum obtained from the data. This allows one 
to calculate the quantities needed to approximate 
the full likelihood function for the binned Ct, using 
the formalism of Bond, Jaffe, and Knox (2000). 

In the second suitable quadratic estima- 

tor of the full sky spectrum in the cut sky variables 
Ct together with it's Fisher matrix is constructed 
via the coupling matrix Mu> and the transfer func- 
tion Fg (Bond, Jaffe, and Knox 1998; Bond et al. 
2001). The underlying power is recovered through 
the iterative convergence of the quadratic estima- 
tor onto the maximum likelihood value as in stan- 
dard Maximum Likelihood (ML) techniques. A 
great simplification and speed-up is obtained due 
to the diagonality of all the quantities involved, 
effectively avoiding the 0(N 3 ) large matrix inver- 
sion problem of the general ML method. The ex- 
tension of the quadratic estimator formalism to 
montecarlo techniques such as MASTER have the 
added advantage that the Fisher matrix character- 
izing the uncertainty in the estimator is recovered 
directly in the iterative solution and does not rely 
on any potentially biased signal+noise simulation 
ensembles. 

The two procedures agree to within a few per- 
cent in the estimated values with the quadratic es- 
timator giving slightly more optimal errors (at the 
5% level) over the sample variance limited range 
of the data. The parameter extraction pipeline 
was run over results from both methods and the 
two were found to agree to within the numerical 
accuracy of the fits. 

The drawback of using naively binned maps in 
the pipeline is that the aggressive time filtering 
completely suppresses the power in the maps be- 
low a critical scale £ c Hivon et al. (2001). This re- 
sults in one or more bands in the power spectrum 
running over modes with no power and which are 
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thus unconstrainable. In pratice we deal with this 
by binning the power so that all the degenerate 
modes lie within the first band 2 < t < 50 and 
regularize the power in the band at DMR power 
in the likelihood analysis. The estimate in the sec- 
ond band 50 > I < 100 will be correlated to this 
regularized value and as such may be also be con- 
sidered to be biased by a prior theoretical input. 
We therefore discard the estimates in the first two 
bands thus avoiding any correlation to the regu- 
larizing scheme used to constrain the power on the 
largest scales. 

An area equivalent to 1.8% of the sky was anal- 
ysed. The region is enclosed in an ellipse with 20 
and 12 degree semi-axes centred at RA = 85 and 
Dec — —46. The data and simulations were pix- 
elised with 7 arcminute pixels (Healpix 7V s id = 
512). The simulations were run at an angular res- 
olution up to £ = 1300. 

Inspection of the BOOMERANG maps shows 
faint stripes of nearly constant declination, (here- 
after "isodec" strips) which vary in amplitude and 
phase between bolometer channels. The striping 
patterns vary on day time scales, and are not re- 
produced in simulated maps made with the same 
scan pattern and best estimated noise correlations 
from the time stream data. 

To eliminate this contaminant, all modes with 
a small k^A (which corresponds to isodec stripes) 
are removed with a Fourier filter. While this 
clearly eliminates isodec stripes, it also filters out 
CMB signal. This is accounted for by applying the 
same filter to the simulated maps in the MASTER 
procedure, so that the effects are included in the 
determination of Fg. The removal of the stripes 
still permits an unbiased estimate of the power 
spectrum of the sky, but does cause a consider- 
able increase in the uncertainties at large angular 
scales. 

The inclusion of several channels is achieved by 
averaging the maps (both from the data, and from 
the Monte-Carlos of each channel) before power 
spectrum estimation. Weighting in the addition 
is by hits per pixel, and by receiver noise at 1Hz. 
Each channel has a slightly different beam size, 
which must be taken into account in the gener- 
ation of the simulated maps. Any inaccuracy in 
assuming a common Bg in the angular power spec- 
trum estimation is then absorbed into Fg. 



The calculation of the full angular power spec- 
trum and covariance matrix for the four good 150 
GHz channels of BOOMERANG (57103 pixels and 
i=a 216, 000, 000 time samples) takes approximately 
1 day running on 8 AMD-athlon 800 MHz work 
stations. 

5. Internal Consistency Tests 

The BOOMERANG observing strategy allows 
for a rich set of internal consistency checks, imple- 
mented as a variety of difference maps in which 
the sky signal should cancel. The power spec- 
tra of these difference maps are sensitive to im- 
proper characterization of the receiver noise, and 
contamination not fixed to the celestial sphere. 
The precision of these difference tests are much 
more powerful than a comparison of the power 
spectra, since the sample variance contribution to 
the power spectrum error bars is proportional to 
the signal found in each bin, which is near zero for 
the difference maps. The results of these tests are 
summarized in Table 2. 

The most powerful of these tests is to difference 
the map made from data acquired while scanning 
at 2dps (the first half of the flight) from data ac- 
quired while scanning at ldps. (the second half 
of the flight). This test is sensitive to solar and 
ground pickup, since between the center of the 
2dps data and the center of the ldps data, the sun 
moves 5° on the sky, and the gondola has moved 
half way around the continent between the time 
centroids of the two maps. The test is also sensi- 
tive to errors in the deconvolution of the transfer 
function of the time-domain signal, since the scan 
speed changes the spatio-temporal mapping of the 
signals. 

This test is performed on each of the five 150 
GHz channels individually. Without filtering out 
the isodec modes as described previously this test 
is failed. With the filtering, 4 of the 5 channels 
pass, and 1 of the 5 channels (150B1) shows a small 
but statistically significant signal. This channel 
is excluded from subsequent analysis. The isodec 
removal is applied to all of the channels included 
in the analysis. 

The 1st half - 2nd half difference test is also 
performed on maps with the four channels com- 
bined. At I < 300, there is a statistically signifi- 
cant residual in the difference map at the level of 
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Table 2 
Internal Consistency Tests. 



Tost 


Reduced % 2 


P> 


B150A ldps - 2dps 


0.91 


0.57 


B150A1 ldps - 2dps 


0.92 


0.56 


B150A2 ldps - 2dps 


1.04 


0.41 


B150B1 ldps - 2dps 


2.73 


7 x 10~ 5 


B150B2 ldps - 2dps 


0.60 


0.91 


4 Ch ldps - 2dps 


1.80 


0.02 


4 Ch Loft - Right 


1.21 


0.24 


(A+Al) - (A2+B2) 


0.61 


0.90 



NOTE. — Reduced x 2 with 19 degrees of freedom 
for internal symmetry tests for BOOMERANG. P> 
gives the probability of obtaining a x 2 larger than 
the one reported. B150B1 fails the test, and is not 
used in the analysis. The '4 Ch' entries combine 
maps from B150A, B150A1, B150A2, and B150B2. 
The ldps-2dps 4 Ch spectrum fails marginally. This 
is dominated by 4 bins centered between / — 150 
and I — 300. The mean signal of these 4 bins is 
50fiK 2 , compared to a signal over the same range of 
ss 4000pif 2 (sec Table 3). 



50/i.ftT 2 . Since the signal only appears in the com- 
bined channel ldps - 2dps analysis, this is consis- 
tent with a noise term which changes slowly on 
the sky and is correlated between channels, such 
as atmosphere. At these angular scales, the CMB 
signal is 4000 [iK 2 . Since the residual signal is 
small compared to the CMB signal, its effects are 
neglected. 

A test for artifacts specific to particular detec- 
tors is made by differencing the map made from 
combining B150A and B150A1 with the combina- 
tion of B150A2 and B150B2, and a test for arti- 
facts due to scan-synchronous baselines is made 
by differencing maps using only the left-going and 
right-going portions of the scans. There is no evi- 
dence of any residual signal in either of these tests, 
which is again consistent with a noise term that 
changes slowly on time scales comparable to the 
scan time. 

6. Foregrounds 

The comparison of the maps at the 4 differ- 
ent frequencies measured by BOOMERANG is a 
powerful tool to test for contamination from fore- 
grounds at 150 GHz. At the resolution frequen- 
cies probed by BOOMERANG, thermal emission 
from interstellar dust grains is expected to be the 



most important foreground (see e.g. Tegmark et 
al. (2000)). Masi et al. (2001) probes the level of 
dust in the BOOMERANG maps by correlating 
BOOMERANG data with a dust template derived 
from the 3000 GHz IRAS/DIRBE maps (Schlegel 
et al. 1999; Finkbeiner et al. 1999) and extrapo- 
lating the dust dominated 410 GHz signal to 150 
GHz using the measured correlations. The de- 
duced power spectrum of dust fluctuations con- 
tributes less than 1% to the power spectrum of 
sky temperature measured at 150 GHz. For this 
reason Galactic dust contamination is neglected in 
the following. 

Radio point sources are another potential form 
of contamination in the maps and angular power 
spectrum. The WOMBAT (1998) radio point 
source extrapolations are used to estimate the ef- 
fects of known radio sources in the BOOMERANG 
fields. The WOMBAT extrapolated fluxes are 
converted to temperature using a Gaussian beam 
that is a good approximation of the beam + pixel 
window function. Assuming that each of the 
WOMBAT sources is in a separate pixel, the rms 
contributed by these point sources to the map is 
calculated. In the Ci power spectrum this should 
show up as a constant Ci = Co contribution, 
which is found by using the effective rms con- 
tributed by a random distribution of point sources, 
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rms 2 = X^(2^+l)C WV(47r), and the known 
beam window function We- For results quoted in 
the units of Table 3, this leads to an estimated 
point source contribution as a function of I of 
160(£/1000)VK 2 . 

However, three quasars are easily identified in 
the maps and removed. The brightest two of these 
are also the two highest flux objects in the Wom- 
bat catalog in our region; the third quasar is the 
eighth brightest in the catalog. Removing only the 
two brightest sources from the catalog and repeat- 
ing the above analytic estimate leads to a predic- 
tion for the contribution of the remaining sources 
of only 85(//1000)VK 2 . 

The power spectrum is evaluated directly from 
the maps before and after removing the three 
quasars. This was done by ignoring pixels within 
0.5° of the quasar position. This cutting induces 
very small additional bin-bin correlations in the 
power spectrum, which are negligible given the 
small area of the cuts. The effect of cutting the 
three quasars is less than 170/xK 2 at all t < 1000. 
This, combined with the analytic estimates above, 
gives us good confidence that the residual radio 
point source contamination is far less than the 
quoted errors at all I. 

7. Power Spectra 

The results are summarized in Table 4 and Ta- 
ble 3, and in Figure 2 and Figure 3. 

The sensitivity of the results to different bin- 
nings is explored. Figure 2 and Table 4 summa- 
rize the results from two independent analysis us- 
ing top-hat binning of width Al = 75, and offset 
and overlapping by 50%. Because of this overlap, 
adgacent points are strongly correlated. This bin- 
ning is not used in the parameter extraction. 

Figure 3 and Table 3 summarized the results 
from an analysis with non-overlapping top-hat 
bins of width Al = 50. The sources of uncertainty 
that are included in the errors quoted in Table 3 
include sample variance and statistical noise. The 
former dominates at £ 600 and the latter at 
higher I. These uncertainties are only weakly cor- 
related. 

Uncertainty in the effective beam-size intro- 
duces an additional uncertainty in the power spec- 
trum that is highly correlated across the spectrum. 
The uncertainty in the effective beam size has con- 




200 400 600 BOO 1000 

I 



Fig. 2. — The angular power spectrum 
of the CMB, as measured at 150 GHz by 
BOOMERANG. The vertical error bars show the 
statistical + sample variance errors on each point. 
There is a common 10% calibration uncertainty in 
temperature, which becomes a 20% uncertainty in 
the units of this plot. The points are also subject 
to an uncertainty in the effective beam width of 
±1.4' (ler). The effect of a la error in the beam 
width would be to move the red points (all up 
together if the beam width has been underesti- 
mated, or all down together if the beam width has 
been overestimated) to the positions shown by the 
black triangles. The blue points would move in 
a similar fashion. The blue and red points show 
the results of two independent analyses using top- 
hat binnings offset and overlapping by 50%. This 
shows the basic result is not dependent on bin- 
ning. While each of the independent spectra (red 
circles or blue squares) are internally nearly un- 
corrected, each red point is highly correlated with 
its blue neighbors, and vice versa. These data are 
listed in Table 4 
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Table 3 
The Angular Spectrum 



Grange ^i±tHc e (fiK 2 ) (ldps-2dps)/2 (fiK 2 ) (left-right) /2 (fiK 2 ) ((A+Al)-(A2+B2))/2 (fiK 2 ) 



( 76 ■ 


- 125) 


3519 


± 


558 


2 


± 


8 


- 11 


± 


11 


4 


± 


4 


(126 


- 175) 


4688 


± 


555 


40 


± 


13 


- 17 


± 


13 


3 


± 


4 


(176 


- 225) 


5592 


± 


548 


23 


± 


14 


- 15 


± 


19 


- 1 


± 


5 


(226 


- 275) 


5699 


± 


486 


45 


± 


19 


- 27 


± 


25 





± 


7 


(276 


- 325) 


3890 


± 


316 


69 


± 


24 


- 44 


± 


33 


- 8 


± 


9 


(326 


- 375) 


2591 


± 


207 


10 


± 


26 


- 32 


± 


44 


- 9 


± 


12 


(376 


- 425) 


1842 


± 


152 


14 


± 


33 


- 98 


± 


55 


- 3 


± 


16 


(426 


- 475) 


2070 


± 


161 


- 58 


± 


37 


- 90 


± 


70 


- 8 


± 


20 


(476 


- 525) 


2267 


± 


174 


24 


± 


53 


- 30 


± 


95 


23 


± 


27 


(526 


- 575) 


2293 


± 


182 


- 9 


± 


68 


- 150 


± 


118 


26 


± 


36 


(576 


- 625) 


2058 


± 


181 


100 


± 


93 


- 161 


± 


155 


- 11 


± 


43 


(626 


- 675) 


1934 


± 


190 


28 


± 


115 


203 


± 


217 


- 23 


± 


56 


(676 


- 725) 


1828 


± 


207 


- 58 


± 


145 


71 


± 


269 


- 32 


± 


71 


(726 


- 775) 


1440 


± 


226 


196 


± 


198 


- 421 


± 


324 


89 


± 


99 


(776 


- 825) 


1920 


± 


288 


-336 


± 


235 


- 808 


± 


411 


160 


± 


131 


(826 


- 875) 


2243 


± 


361 


-211 


± 


317 


- 73 


± 


580 


176 


± 


171 


(876 


- 925) 


1752 


± 


428 


- 94 


± 


437 


- 613 


± 


757 


- 23 


± 


217 


(926 


- 975) 


985 


-j- 


506 


- 78 


-j- 


591 


- 607 


± 


1013 


-458 


± 


278 


(976 ■ 


■ 1025) 


502 




627 


-128 




800 


-1370 


± 


1347 


- 82 


± 


395 



Note. — The spectrum of the CMB, as used in the parameter extraction. The spectrum is further subject to an 
overall 10% calibration uncertainty, and a 1.4' effective beam uncertainty. The spectrum of the all-channel consistency 
tests are also given. Adjacent bins are weakly correlated 



tributions from uncertainty in the physical beam 
and from uncertainty in the rms amplitude of the 
pointing jitter. These combine to produce an un- 
certainty in the effective beam of ±13%. This un- 
certainty is not included in the errors quoted in 
Table 3, as its effect is to produce an overall tilt 
to the spectrum. The amplitude of the tilt cor- 
responding to the 1 sigma uncertainty that is as- 
signed to the effective beam width is illustrated 
in Figure 2. This uncertainty is included in the 
parameter estimation outlined in the next section. 

Uncertainty due to instrumental and atmo- 
spheric artifacts in the maps are small, based on 
the internal consistency tests, and are neglected. 
Similarly, contamination of the maps by both dif- 
fuse and compact astrophysical foregrounds are 
also negligible with respect to the other uncertain- 
ties and are neglected. 

8. Cosmological Parameters 

The angular power spectrum shown in Figure 3 
can be used in conjunction with other information 
to determine cosmological parameters. In this pa- 
per, the parameter extraction methods and tools 
described in detail in Langc et al. (2001) are used. 



Specifically, the relative agreement between 
these data and theoretical predictions over a broad 
7-dimensional cosmological parameter space is ex- 
plored. This parameter space is appropriate for 
models with adiabatic initial conditions (e.g. in- 
flationary models). The resolution of the second 
peak virtually rules out alternative models such 
as defect based scenarios which predict a single 
broader peak with no secondaries (Turok et al. 
1998; Contaldi et al. 1999). 

Parameters explored include those describing 
energy densities, including the total energy den- 
sity f^tot, the vacuum energy density f^A, and the 
physical densities of baryons and cold dark mat- 
ter, flbh 2 and fl c h 2 respectively. The power spec- 
trum of initial adiabatic density fluctuations is de- 
scribed by a normalization In C\o and a power law 
exponent n s . The effects of recent reionization of 
the universe on the observed angular power spec- 
trum, parameterized by the optical depth to de- 
coupling, tc, are also explored. For this parame- 
ter, the liklihood does not fall sufficiently by the 
edge of the explored parameter range (r c < 0.5) 
to produce significant limits, though a preference 
for low values of tau c are evident in Figure 4. 

Given the data, likelihoods as a function of 
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Table 4 
The Angular Spectrum 



Grange ^lc e (piK 2 ) 



( 50 - 


113) 


3035 


± 


557 


( 75 - 


150) 


3776 


± 


428 


(112 - 


187) 


4481 


± 


416 


t^iou - 


225) 


5380 


± 


429 


(187 - 


262) 


5810 


± 


413 


(225 - 


300) 


5245 


± 


345 


(262 - 


337) 


4056 


± 


257 


(300 - 


375) 


2942 


± 


184 


(337 - 


412) 


2218 


± 


140 


(375 - 


450) 


1861 


± 


119 


(412 - 


487) 


1992 


± 


123 


(450 - 


525) 


2424 


± 


138 


(487 - 


562) 


2443 


± 


142 


(530 - 


605) 


2520 


± 


162 


(567 - 


642) 


2298 


± 


160 


(600 - 


675) 


1868 


± 


144 


(637 - 


712) 


1858 


± 


154 


(675 - 


750) 


1696 


± 


163 


(712 - 


787) 


1560 


± 


179 


(750 - 


825) 


1736 


± 


211 


(787 - 


862) 


2021 


± 


250 


(825 - 


900) 


2172 


± 


292 


(862 - 


937) 


1847 


± 


333 


(900 - 


975) 


1174 


± 


377 


(937 - 


1012) 


762 


± 


437 


(963 - 


1038) 


499 


± 


623 


NOTE.- 


—The spectrum 


of the 


CMB, as 


shown in 


Figure 2 


. These 


are the 


the results of 


two inde- 


pendent 


analyses 


using 


Al 


= 75 


top-hat 


binnings 


offset 


by 


Al = 



75/2 and thus overlapping by 50%. 
While each of the independent 
spectra are internally nearly uncor- 
rected, each red point is highly 
correlated with its neighbors from 
the other binning. This binning is 
not used in the parameter extrac- 
tion. Rather, the non-overlapping 
(and thus only weakly correlated) 
Al — 50 top-hat binning listed in 
Table 3 is. The spectrum is further 
subject to an overall 10% calibra- 
tion uncertainty, and a 1.4' effective 
beam uncertainty. 
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Table 5 

Results of Parameter Extraction 



Priors 


fltot n s flbh 2 Qcdm 


h 2 


n A 




ft 


Ago 


Weak only 


i no0.06 n ncO.lO n n22 0.004 n 10O 
1 - UZ 0.06 U ' yD 0.09 U ' UZZ 0.003 U.IOq 


05 
05 


Co.5i8:lg) 


(0.5l8:lg) 0.07° : °J 


(0.56°:1°) 


i5.2i:» 



LSS 
SNla 

LSS & SNla 
ft = 0.71 ± 0.08 



1.02" 
1.02° 
0.99" 
0.98°, 



0.97^ 
0.99°, 
1.03" 
0.972 



029"""' 1 
023 004 

U.U^Og QQ 4 

U.UZOq 003 

022 004 



0.13! 
0.10, 
0.14 
0.14' 



o 02i luula 

0.021°-»f 
0.022° 
0.023!! 



11.(14 

0.04 
0.01 
0.01 
0.01 

81 



0.5511 
0.73° 
0.65o 
0.62°, 



u -* y 0.12 

u - ol 0.06 
34°°* 

40°- 13 

(0.48^2) 

0.38° : °? 

0-33 : 1 

0.33° : °i 



fl7 u u ^ 
U.U/g.02 

06° ° 3 

n 05°° 2 

U.UOn.Q2 
°.° 5 0.02 
U - UD 0.02 

o.o5 1; 

0.05 If 

o.o5°°; 



0.56 (,n 

61°° 9 
67°° 9 
(0-65g : g|) 
(0-61-il) 



i5.oi:| 

15 9 2 5 

13.71J 

13 7 16 
' 1.6 
oil. 6 
•'0.6 
5 0.7 



Flat 

Flat & LSS 
Flat & SNla 
Flat, LSS & SNla 



(1.00) 
(1.00) 
(1.00) 
(1.00) 



0.95 
0.98 
0.98 
1.03! 



0.08 
0.10 
0.07 
0.11 



111!!! 



J 0.003 



0.13! 
0.13, 
0.12, 
0.13 



0:01 
0.01 



(0.57°:^) 
0.62°:°* 
0.68°° 4 
0.66°, 



u 0.06 



0.62 
0.66 
0.66 



00 



14-3 ' 
14.5| 
14.0° 
14.0° 



Note. — Results of parameter extraction using successively more restrictive priors, following Lange et al. (2001). The confidence 
intervals are la. The quoted values arc reported after marginalizing over all other parameters. For the primary database parameters, 
16% and 84% integrals are reported as ±1(T errors. For O m , O^,, h, and Age, which are functions of the other parameters, the 
mean and standard deviation over the distribution arc reported. All entries arc subject to a weak prior in which only models 
with 0.45 < h < 0.90 and age > 10 Gyr are considered. The LSS (Bond & Jaffe 1999) and SNla supernovae (Riess et al. 1998; 
Pcrlmutter et al. 1999) priors arc as described in Lange et al. (2001). The strong h prior is a Gaussian with the stated 1<7 error. 
Parentheses are used to indicate parameters that did not shift more than l-a or have their errors reduced by a factor of two upon 
the inclusion of the CMB data, compared with an analysis using the priors only. Thus, in these cases the parameter range reflects 
the choice of prior, rather than a constraint by the CMB. The age column is in units of Gyr. 



theoretically predicted power spectra are calcu- 
lated throughout this parameter space. For ev- 
ery comsological model, beam and calibration un- 
certainties add two additional parameters. We 
approximate the possible correction to the beam 
of effective width u> by a gaussian form Wg(u> + 
5uj)/W t {uj) = e -^+i)^ ; with 5lo assumed to be 
gaussian distributed with the standard deviation 
corresponding to 1.4'. Effectively, every theoreti- 
cal spectrum is multipled by e 2 ^+ 1 ) w < 5 "\ Calibra- 
tion uncertainty of 10% adds to the variance of 
lnCio- We do not calculate likelihoods on a grid 
in beam- width - calibration space, rather for every 
model we search for a maximum of likelihood in 
5w and lnCio, calculate the curvature of the like- 
lihood near this maximum, and marginalize over 
the beam and the amplitude by integrating likeli- 
hood in the gaussian approximation. The best fit 
value of lnCio an d it's variance is used when com- 
bining Boomerang predictions with other data. 

Parameters are constrained individually by 
marginalizing over all others, including two de- 
scribing the calibration and beam window func- 
tion uncertainties. The 16% and 84% integrals 
are reported as ±lcr errors. Other quantities such 
as the Hubble constant and the age of the Uni- 
verse are derived from those used to define the 
parameterization, using the mean and variance 



over the posterior distribution. The details of the 
discrete numerical database used for this process, 
including the limits and values of each parameter 
tested, and exact prescription used for calculat- 
ing likelihoods and extracting confidence limits, is 
fully described in Lange et al. (2001). 

Before marginalization, the calculated likeli- 
hood for each model is multiplied by the likelihood 
derived from a series of "prior probabilities", or 
priors, which represent knowledge from other cos- 
mological measurements. All results considered 
here have a "weak h + age" top-hat prior ap- 
plied (hereafter simply the "weak prior") which 
eliminates models where the universe is younger 
than 10 Gyr, and limits the Hubble constant, 
H = 100 h km s^Mpc -1 , to 0.45 < h < 0.9. 

Applying stronger priors in conjunction with 
this weak prior exercises the ability of the CMB 
data to combine with other measurements (or 
theoretical prejudice) to significantly narrow the 
parameter confidence intervals. Considered here 
is the impact of applying a stronger constraint 
on h, constraints derived from measurements of 
large scale structure (LSS) (Bond & Jaffe 1999), 
results from recent measurements of type la super- 
novae (Riess et al. 1998; Perlmutter et al. 1999), 
and the theoretical bias that fitot = 1. 

The parameter estimates, given these various 
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Fig. 3. Selected best fit models normalized to 

the best compromise amplitude between COBE-DMR and 
BOOMERANG are shown overlayed on the BOOMERANG 
spectrum. The upper panel shows the points plotted as listed 
in table 3. The best-fit models for the "weak" and "no pri- 
ors" cases coincide (blue, solid curve) with Qtot — 1.05, Hq = 
50, fi A = 0.5, wj = 0.020, w c = 0.120, r e = 0, n s = 0.925, t = 
15.8Gyrs. Strong Hubble prior gives the best fit model with 
parameters Qtot — 1-, Ho — 68, Q\ — 0.7, Ub — 0.020, u c — 
0.120, r c = 0,n 3 = 0.925, to = 13.8Gyrs. The full analy- 
sis takes into account the calibration and beam uncertainties 
which best fit models take advantage of. This explains the 
apparent mismatch between some of the models in the up- 
per panel and the plotted central values of Boomerang band 
powers. The green (dash-dot) curve is the best fit model 
(O tot = 1.15, Ho = 42, Q A = 0.7, u b = 0.020, w c = 0.060, t c = 
0.2, n s — 0.925, to — 20Gyrs) when both beam and calibra- 
tion uncertainties are switched off. The model fits closely the 
central values as expected. To demonstrate the effect of beam 
and calibration unccrtanties, in the lower panel the data points 
have been replotted with a 4% decrease in calibration (0.4er) 
and a 0.5 arcminutc change in beam size (0.4<r). The plot 
makes it clear that the best-fit conventional CDM models arc 
indeed good fits to the data, once these uncertainties are cor- 
rectly accounted for. 



combinations of priors, are shown in Table 5, with 
marginalized likelihood curves for several impor- 
tant parameters given in Figure 4. 

Lange et al. (2001) reports a family of models 
within the database that provide good fits to the 
angular power spectrum up through I ~ 600, but 
represent very young, high baryon density, very 
closed models outside the normal realm of con- 
sideration in modern cosmology. In Lange et al. 
(2001), the weak prior was used to keep these mod- 
els from affecting the parameter estimates. Using 
the power spectrum shown in Figure 3, but limit- 
ing consideration to the points with I < 600, sim- 
ilar behavior is exhibited; the very young, high 
baryon density, very closed models dominate the 
fits, but can be eliminated by the weak prior. As 
in Jaffe et al. (2001), adding the higher I points, 
which exclude models with either a very high or 
very small third peak, eliminates the high baryon 
density models. It does not eliminate the closed 
models in the absence of the weak priors. Thus, 
just as in Lange et al. (2001), cases with the weak 
priors applied arc the focus of the extraction. 




0.05 2 4 0.6 0.6 0.8 1 1.2 1.4 

n b h 2 n c h z n . 



Fig. 4. — Likelihood functions for a subset of the 
priors used in Table 5. flbh 2 , and n s are well con- 
strained, even under the "whole database" case, 
and are insensitive to additional priors, fitot is 
poorly constrained over the whole database, but 
when the weak priors are applied, it becomes sta- 
bly consistent with the flat case. With the weak 
priors, il\ and Q c h 2 are poorly constrained, but 
become significant detections with the addition of 
the other priors considered. For all cases, r c is 
poorly constrained, but does prefer low values. 

As is apparent in Figure 4 and Table 5, il^h 2 
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and n s are well localized for all choices of priors. 
The range for flbh 2 is very consistent with deter- 
minations based on Big Bang Nucleosynthesis and 
measurements of light element abundances (Tytler 
et al. 2000), while n s is localized near unity, consis- 
tent with inflation-based models. Similarly, once 
the weak priors arc applied, f2 to t is well con- 
strained and consistent with a flat universe. f^A 
and n c h 2 are poorly constrained for the weak prior 
case, tc favors low values, but not at a statisti- 
cally significant level. 

These limits agree well with those found in 
Lange et al. (2001), with one exception. There, 
the lcr range for Vt^h 2 (with the weak prior) is 
0.036 ± 0.005. Considering only the points in 
Figure 3 with t < 600 results in a lcr limit of 
0.027 ±0.005. This shift is presumably due to the 
improvement of the pointing solution and smaller 
error bars compared with BOO. The table reflects 
the estimate made using all points up to I = 1000; 
the addition of the information contained in the 
high-^ points has moved the confidence intervals 
still further. 

Having found tight limits on flbh 2 , fi to t ; and n s 
with the weak priors, the impact of adding other 
priors can be discussed, to see whether these re- 
sults are stable (i.e., consistent with the prior) and 
whether any new parameters can be localized. The 
table indicates that, in fact, these three parame- 
ters are very stable to the addition of the LSS, 
SNla, strong h and the theoretically motivated 
ft tot = 1 priors, in various combinations. This 
insensitivity to choice of priors is a powerful indi- 
cation that the cosmology is consistent and that 
these parameters have been robustly determined. 

The table also shows that the CMB data can be 
combined with these priors independently to make 
statistically significant determinations of ftA and 
ft c /i 2 . While the confidence intervals shift some- 
what depending on the chosen prior, the rough 
agreement among these three results, giving ftA ~ 
0.65 and Q c h 2 ~ 0.12, is very compelling. These 
LSS and SNla results are similar to those found 
in Lange et al. (2001). For the first time, how- 
ever, the combination of CMB data with just a 
strong limit on h is powerful enough to yield such 
a detection. 

The table also shows a consistent story for the 
age of the universe and for the Hubble Constant. 
For the prior combinations where the limits on h 



are dominated by the CMB data rather than the 
priors, the extraction finds h ~ 0.65, with ages 
between 13 and 15 Gyr. These quantities are most 
strongly constrained by the CMB data along with 
the combination of the SNla prior and flatness. 

9. Conclusions 

A measurement of the angular power spectrum 
of the CMB, characterized by a series of harmonic 
peaks, has been presented, confirming the exis- 
tence of this unique signature of acoustic oscil- 
lations in the early universe. This is an impor- 
tant confirmation of standard adiabatic models 
of structure formation, and thus of the process 
of constraining cosmological parameters based on 
increasingly precise measurements of the position 
and amplitude of these peaks by current and fu- 
ture CMB experiments. 

The precision and extent of the angular power 
spectrum that is reported here already allow an 
accurate determination of several cosmological pa- 
rameters with the assumption of only weak astro- 
physical priors. Assuming 0.45 < h < 0.90, the 
CMB data tightly constrain the values of fttot and 
n s to lie close to unity and tightly constrain ft;, to 
a value consistent with BBN. 

Adding constraints from observations of Large 
Scale Structure and of type la supernovae, yields 
a value for the Hubble constant of h = 0.67±0.09, 
that is in good agreement with the HST key 
project final value of 0.72 ± 0.08 (Freedman et al 
2000). These data also provide compelling evi- 
dence for the existence of both dark matter and 
dark energy. Including LSS, SNla or a prior on 
the Hubble constant of h = 0.71 ± 0.08 each yields 
ft A ~ 2/3 and ft m w 1/3. 
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port from Kathy Deniston, and superb field and 
flight support from NSBF and the USAP person- 
nel in McMurdo. 
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